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INTRODUCTION 

It  was  first  reported  by  Stanyukovich  (Ref.  1)  that  the 
launching  of  a  projectile  from  a  gas  gun  can  be  accomplished 
with  constant  pressure  at  the  base  of  the  projectile  if  a 
particular  kind  of  propellant  flow  is  utilized  in  the  gun 
barrel.  Stanyukovich  derived  the  equation  of  the  required 
flow  from  the  assumption  that  the  flow  velocity  is  a  function 
of  time  only.  Independently  of  Stanyukovich,  Curtis  (Ref.  2) 
and  the  author  (Ref.  3)  arrived  at  the  same  analytical  solu¬ 
tion,  but  used  different  approaches.  Curtis  derived  the 
solution  from  an  analogy  with  the  earth's  atmosphere,  the 
author  from  the  necessary  and  sufficient  condition  to  main¬ 
tain  constant  pressure  of  each  element  of  a  fluid  in  motion. 
Various  Investigators  have  elaborated  on  the  solution  and 
suggested  modifications  of  two-stage  gas  guns  to  approximate 
the  conditions  for  a  constant  base  pressure  launching  of  a 
projectile  (Refs.  2  through  6). 

In  the  present  report  the  parti 'T.»lar  flow  required  in 
the  gun  barrel  will  be  studied  by  the  method  of  characteristics. 
The  investigation  will  not  be  restricted  to  the  propellant 
flow  in  the  gun  barrel,  but  will  be  extended  to  the  more  gen¬ 
eral  flow  phenomenon,  which  has  been  termed  "constant  accel¬ 
eration  flow"  by  the  author.  To  facilitate  the  presentation 
by  the  method  of  characteristics  the  production  of  a  constant 
acceleration  flow  by  a  simple  means  will  be  investigated  in 
more  detail.  The  constant  acceleration  flow  equations  will 
be  reviewed  briefly,  for  reason  of  reference,  and  some  relations 
will  be  derived  which  have  not  been  given  elsewhere. 
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CONSTANT  ACCELERATION  FLOW  EQUATIONS 


Assualng  a  one  diaensional  flow  in  the  absence  of  vis¬ 
cosity  and  heat  conduction,  the  author  has  shown,  in  reference 
3,  that  the  pressure  of  each  fluid  elesent  can  only  reaain 
constant  in  tiae  if  no  relative  velocity  variation  exists 
throughout  the  flow.  This  is  obviously  the  case  for  a  fluid 
at  rest  or  in  constant  velocity  aotion.  But  the  fluid  also 
exhibits  no  relative  velocity  variation  if  the  saae  force  is 
acting  on  the  unit  aass  throughout  the  fluid  at  each  instant 
of  tiae.  In  the  absence  of  body  forces,  pressure  tlaes  area 
forces  are  active  only  and  it  follows  that  the  force  is  con¬ 
stant  and  thus  produces  a  flow  of  constant  acceleration. 

The  proof,  that  the  "no  relative  velocity  variation"  is 
a  necessary  and  sufficient  condition  for  the  pressure  of  each 
fluid  eleaent  to  reaain  constant  in  tiae,  follows  froa  the 
continuity  equation 


*  p 


Du 

Dx 


O 


(1) 


with  du/dxm  o  one  obtains  and  it  can  be  shown 

(see  Ref.  3)  that  G^P>/utmO  **  *  consequence  of  c^P/c^tmO 
Independent  of  both  the  process  by  which  the  flow  has  been 
generated  and  the  equation  of  state  of  the  fluid  used.  For 
the  strictly  one-dlaensional  flow,  Euler's  equation  of  notion 
becoaes  with  Du/$x  *  0 


2u 

JT 


0( 


(2) 


where  oc  is  a  constant. 

The  constant  acceleration  flow  exhibits  a  pressure  gra¬ 
dient  which  does  not  vary  with  tiae  and  a  constant  temporal 
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velocity  gradient.  A  temperature  gradient  is  present  in  the 
flow  unless  the  flow  has  been  produced  by  an  isothermal  pro¬ 
cess.  Using  Lagranglan  coordinates  one  can  easily  show  that 
none  of  the  properties  of  each  fluid  element  changes  with  time. 
The  integration  of  the  equation  of  motion  leads  to 

(XX  +  ^  f ft J 

where  pQ  is  a  reference  pressure  and  f(t)  is  the  time  dependent 
constant  of  integration.  Since  neither  p  nor  p  depend  on  time, 
differentiation  with  respect  to  time  gives 

f'(t)  =  *  -occL-  cx*£ 


if  u  -  o  for  t  -  o.  The  general  equation  of  a  constant  accel¬ 
eration  flow  then  becomes 


f*  *  **  -  O 

P 


(3) 


For  an  isentropic  p-p  relation  of  an  ideal  gas  one  obtains 
from  Eq.  3 

(4) 


x  r-  / 

/>/,,  Jr  ,  r-/  /». 

(  /P°J  = 7  t*: 


OCX  + 


-2<r  /»• 


where  pQ  and  pQ  denote  their  values  at  x  -  o  and  t  -  o.  Intro¬ 
ducing  the  speed  of  sound  in  Eq.  4  we  have 

a  "*«•  **  ^  (5) 

where  a  -  aQ  for  x  -  o  and  t  -  o.  This  equation  has  first 
been  derived  by  Stanyukovlch  (Ref.  1) 
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CONSTANT  ACCELERATION  FLOW  PROPERTIES 


The  simplest  and  most  illustrative  way  to  produce  a  con¬ 
stant  acceleration  flow  is  to  set  a  horizontal  pipe  enclosing 
a  certain  amount  of  fluid  into  motion  of  constant  acceleration. 
Depending  on  the  degree  of  acceleration  and  the  length  of  the 
pipe,  the  transition  of  the  enclosed  gas  from  rest  into  the 
accelerated  motion  may  or  may  not  lead  to  the  formation  of  a 
shock  wave.  In  order  to  avoid  this  it  is  assumed  here  that 
the  acceleration  <x  is  approached  sufficiently  slowly.  The 
adjustment  of  the  enclosed  gas  to  constant  acceleration  motion 
then  has  taken  place  isentropically.  To  describe  the  condi¬ 
tions  within  the  gas  in  Eulerian  coordinates,  the  time  will  be 
counted  from  the  instant  the  acceleration  oc  is  reached.  The 
origin  of  the  coordinate  system  is  placed  at  the  location  of 
the  reference  pressure  p0  (Eq.  4)  at  time  t  -  o  and  held  sta¬ 
tionary.  The  positive  x-  axis  points  into  the  direction  of 
motion  of  the  fluid.  The  conditions  in  each  cross  section  of 
fluid  are  then  given  as  functions  of  position  and  time  by  Eqs. 

4  and  5.  Simpler  expressions  of  the  constant  acceleration 
flow  are  obtained  by  means  of  Larrangian  coordinates.  The 
coordinate  system  then  moves  with  the  fluid  and  the  origin  is 
attached  to  the  cross  section  which  exhibits  the  reference 
pressure  pQ  after  the  acceleration  oc  is  reached.  Since  none  of 
the  gas  properties  change  with  time  the  conditions  throughout 
the  fluid  are,  according  to  Eqs.  4  and  5,  given  by 


and 


/  ~ 


Po_ 

<r  jo0 


cx  S 


(6) 


c*0z  - 


(7) 
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The  coordinate  s,  used  here  for  reason  of  distinction,  repre¬ 
sents  the  distance  of  a  given  fluid  particle  from  the  reference 
cross  section  at  the  time  cx  is  reached.  The  values  of  s 
increase  in  the  direction  of  the  accelerated  notion.  In 
describing  the  flow  phenonenon  further,  we  will  use  Lagrangian 
presentation  for  reason  of  convenience. 

The  above  equations  indicate  that  the  pressure  p,  the 
density  p,  and  the  sound  speed  become  zero  for  a  certain 
coordinate  a*.  This  cross  section  is  given  by 

♦  jr  __  ct0 

3  *  7^”  Pa  “  7 

A  vacuun  thus  exists  in  that  part  of  the  pipe  for  which  s  ■  s*. 


The  temperature  distribution  is  obtained  from  Eq.  7  as 


T- 


(9) 


where  cp  is  the  specific  heat  at  constant  pressure.  The  tem¬ 
perature  thus  decreases  linearly  with  s  and  becomes  zero  for 
s*  m  c^To/^  ,  The  enthalpy  also  decreases  linearly  with  s 
according  to 


-  &  o  ~  3 


(10) 


with  zero  enthalpy  for  the  particle  at  s*  -  h Q/ (X . 

The  variation  cf  pressure,  density  temperature,  and  sound 
speed  as  a  function  of  s  are  shown  in  Fig.  1.  The  pressure 
and  density  approach  gradually  with  horizontal  tangent  the 
value  zero.  The  speed  of  sound  according  to  its  parabolic 
form  of  variation 
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('r-O  *  SJ 


(id 


decreases  abruptly  with  vertical  tangent  to  zero. 

For  values  of  s  >  s*  the  enthalpy  and  temperature  becomes 
negative,  the  pressure,  density,  and  sound  speed  become  imagi¬ 
nary.  A  constant  acceleration  flow  cannot  physically  be  real¬ 
ized  beyond  s*.  It  will  be  shown  later  that  the  zero  density 
particle  path  plays  the  particular  role  of  an  envelope  of 
characteristics. 

The  relation  will  now  be  derived  between  the  conditions 
in  the  pipe  before  and  after  the  acceleration  is  reached. 

This  could  be  done  by  evaluating  the  average  density  p  from  an 
integration  of  Eq.  6.  The  average  density  can,  however,  be 
derived  with  less  effort  by  direct  integration  of  Euler's 
equation  at  constant  time.  It  follows  then  for  a  pipe  of 
length  s  that 

/£>„-/£>  =  tx  /  pat*  -  -  °rV°-y 

where  mff  is  the  mass  of  the  enclosed  gas.  Eq.  12  is 
independent  of  the  p-p  relation.  The  above  equation 
particularly  simple  if  the  pipe  is  of  length  s+.  It 
then  that 

P„  -  <xp  -s  * 


valid 

becomes 

follows 

(13) 


and  one  obtains  by  replacing  s*  by  its  value 


r 

r- ' 


p 


(14) 
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It  is  now  interesting  to  determine  at  what  distance  s  the 
average  density  p  occurs.  At  this  point  all  state  variables 
will  be  the  same  as  they  had  been  before  the  pipe  was  acceler¬ 
ated.  Eq.  14  and  Eq.  6  give 


From  Eq.  14 ,  15,  9  and  6  one  obtains 

04  6  S* 

o.  *op0 

O.  S4  7Z 
0  .2  2 


For  %  Par-  y  .  % 

s  -  O.  39  s*  s  = 


p 

=  O  29  Po 

P 

T 

*  0.6/71 

f 

P 

-  0./7  po 

P 

(15) 


To  illustrate  the  generation  of  a  constant  acceleration 
flow  further,  it  will  be  assumed  that  the  flow  is  produced  by 
erecting  a  pipe  enclosing  a  certain  amount  of  gas  from  a  hori¬ 
zontal  to  a  vertical  position.  The  acceleration  then  equals 
the  acceleration  of  gravity  g.  A  Lagrangian  coordinate  system 
will  be  placed  with  the  origin  at  the  bottom  of  the  pipe  and 
the  positive  s-  axis  vertically  up.  Assume  the  pipe  has  been 
filled  in  a  horizontal  position  with  air  of  one  amagat  at  a 
temperature  of  T  -300°K.  One  finds  then  that  the  pressure 
becomes  zero  at  the  upper  end  of  the  pipe  if  its  length  is 
s*  -  31km.  The  density  at  the  bottom  enclosure  of  the  pipe 
has  increased  to  3.5  amagat,  the  temperature  to  495°K,  and 
the  pressure  to  5.8  atmospheres.  The  distance  s,  where  the 
density  of  one  amagat  would  be  found,  is  at  a  height  of  s-12km. 
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This  art! f leal  atmosphere,  which  was  created  by  an  lsen- 
tropic  adjustment  from  the  horizontal  to  the  vertical  position 
of  the  pipe,  differs  noticeably  from  an  ordinary  atmosphere. 

If  by  heat  conduction  the  same  temperature  could  be  established 
throughout  the  enclosed  gas,  the  conditions  then  become  similar 
to  those  of  the  surrounding  atmosphere.  The  "flow"  now  present 
in  the  pipe  represents  an  isothermal  constant  acceleration  flow. 
The  existence  of  this  flow  in  the  earth's  atmosphere  was  proved 
by  the  author  (Ref.  3)  employing  Eulerian  coordinates.  The 
properties  of  an  Isothermal  constant  acceleration  flow  will  now 
be  given. 

With  p/p  •  const,  the  solution  of  the  general  equation  of 
a  constant  acceleration  flow  (Eq.  3)  becomes  In  Eulerian 
coordinates 

*•//>.  ^  +  «x  -  -  O  (16) 

In  Lagrangian  coordinates  one  obtains 

<17) 

The  Lagrangian  presentation  shows  bast  the  similarity  to  an 
isothermal  atmosphere  where  the  pressure  decreases  exponen¬ 
tially  with  increasing  height.  The  significant  difference 
of  the  constant  acceleration  flows  involving  either  an  isen- 
tropic  or  an  isothermal  process  is  that  in  the  isothermal 
case  the  pressure  becomes  zero  only  at  infinity. 
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CONSTANT  ACCELERATION  FLOW  CHARACTERISTICS 

The  equations  of  the  characteristics  can  easily  be  derived 
since  a  and  u  are  known  as  functions  of  x  and  t.  We  write  the 
characteristics  in  the  fora  -t£-/  U  +  a  -  k  .  The  families 
of  curves  which  present  the  net  of  the  characteristics  are 
then  generated  by  using  k  as  *  parameter.  With  u  -  o  for  t  -  o 
it  is  most  convenient  to  use  as  a  parameter  ax  which  represents 
the  local  speed  of  sound  along  the  x-  axis.  The  value  of  ax 
is  given  by  Eq.  5.  We  have 

a,  .  a„  {/  -  -^4  *  *)  (18) 

where  aQ  equals  the  sound  speed  of  the  fluid  element  initially 
at  position  x  -  o.  The  equations  expressing  the  net  of  char¬ 
acteristics  become 


and  with 


U  ±  a 


(19) 


a 

u 


v-' 


cx  £ 


one  obtains  by  rearrangement 


*  \  \  4 

~  (3- fjoc 


fs-yjcr*-  2<?xxj 


(20) 


The  characteristics  are  parabolas.  The  nxes  of  the  parabolas 
are  all  parallel  to  the  x-  axis.  With  x L  7  and  ^  u  (  ) 
the  parabolas  become 


{ $ -r) 


(21) 


All  parabolas  are  of  this  identical  shape. 
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Eliminating  ax  between  the  expression  for  x*  and  t'  we 
derive  an  equation  which  describes  the  locus  of  the  vertices 
of  the  parabolas  as 

*0r-)  /_5»«  _  r) 

C*f3-yj  ( /  (22) 

Since  the  vertices  of  the  parabolas  represent  for  the  charac¬ 
teristics  the  points  of  zero  slope,  the  equation  (22)  indicates 
where  sonic  conditions  occur  in  the  flow.  The  above  parabola 
will  therefore  be  called  the  sonic  parabola.  Its  equation 
could  also  have  been  derived  from  the  condition  u+ar*0. 

The  sonic  parabola  Intersects  the  x-  axis  at  a  distance 
* 

of  x  sr  - - f-—  where  ax  -  o.  The  intersection  with  the 

Cr-'Jo<  - 

t  -  axis  is  at  t  -  t  .  At  this  time  sonic  conditions 

occur  in  a  constant  acceleration  gun  at  the  gun  barrel  entrance. 
The  particle  path  describes  a  parabola  of  the  form 

(23) 

where  xQ  indicates  the  initial  position  of  the  particle. 

Two  particle  paths,  of  initial  position  x  -  o  and  x  - 

,  the  sonic  parabola  and  two  characteristics  are 
shown  in  Fig.  2.  A  schematic  net  of  characteristics  is  illus¬ 
trated  by  Fig.  3.  Notations  have  been  omitted  in  this  figure 
so  as  not  to  distract  from  its  Op-Art  value. 

We  observe  that  the  system  of  characteristics  is  symmetri¬ 
cal  to  the  x-  axis  and  that  the  separation  of  neighboring 
characteristics  reduces  with  increasing  x  value.  The  inter¬ 
section  points  of  infinitely  closely  spaced  characteristics 
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define  an  envelope  with  the  equation 

t*  -  “*/#/’*“  (24) 

This  envelope  is  a  parabola  whose  vertex  is  at  x  -  /^-/ Jcx. 

where  ax  *0.  Its  equation  describes  the  path  of  the  particle 
initially  at  a  location  where  ax  -  o.  The  gas  element  which 
moves  along  this  path  exhibits,  as  has  been  shown,  zero  pres¬ 
sure,  temperature,  density,  and  sound  velocity.  The  physically 
realizable  segment  of  the  x,  t  plane  in  which  a  constant  accel¬ 
eration  flow  produced  by  an  isentropic  process  can  exist  is 
thus  the  above  zero  density  particle  path  and  the  x-  axis.  To 
the  right  of  the  zero  density  particle  path  the  sound  velocity 
assumes  Imaginary  values.  The  separation  of  neighboring  char¬ 
acteristics  increases  again  in  this  region. 

For  the  particular  case  of  a  constant  acceleration  gun 
with  projectile  position  at  x  -  o  for  t  -  o,  the  constant 
acceleration  flow  exists  only  to  the  left  of  the  particle  path 
originating  at  x  -  o.  Sound  velocity  is  reached  at  the  base 
of  the  projectile  where  the  path  of  the  particle  originally  at 
x  -  o  intersects  the  sonic  parabola.  It  occurs  at  the  time 
t  m  / <x  .  If  x  -  o  is  the  position  of  the  barrel  entrance, 
sonic  conditions  occur  there  at  the  time  indicated  by  the  inter¬ 
section  of  the  sonic  parabola  with  the  positive  t-  axis.  This 
time  is  t  •  / /4-<r  .  Th^  propellant  flow  will  be  restricted 

to  u  ■  a  from  now  on  for  a  chambered  gun. 

The  outstanding  feature  of  the  constant  acceleration 
flow  is  especially  well  illustrated  by  presenting  the  charac¬ 
teristics  in  Lagrangian  coordinates.  We  select  again  ax  as 
the  parameter.  The  equations  of  the  characteristics  for  the 
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isentropicmlly  produced  constant  acceleration  flow  follow  from 
the  relations  which  express  u  and  a  as  functions  of  the  particle 
coordinate  s  and  time  t.  We  have  u  -  t  and  a(s,t)~ 
ao  { / -  c*sj/x  since  in  Lagrangian  coordinates  -O  . 

We  obtain  then 


-2. 


(25) 


or  after  rearrangement 

-  )*  1  -s7-  fs+-s7 

+  fr-/JcxS  fr-'J*  J  fr-'J*  *-  V  (26) 

etc9' 

where  a*  has  been  replaced  by  a8  and  ^ ^  by  s*.  Setting 
f  J  ■  t*  and  L  J  “  s*  the  parabolas  can  be  transformed 
into  one  parabola  with  the  equation 


Zf 


'a 


* 

rr-')<x 


(27) 


The  vertices  form  the  envelope 


3^ 


C*o 


X 


* j  o< 


(28) 


For  this  s  value  we  have  a8f  pressure,  density,  and 
temperature  equal  to  zero.  To  the  right  of  the  envelope  the 
sound  velocity  assumes  Imaginary  values.  The  positive  sign  in 
equation  (25)  is  satisfied  by  the  lower  branch  of  the  parabola 
of  equation  (26)  and  its  extension  as  the  upper  branch  into 
the  region  of  imaginary  sound  speeds,  the  negative  sign  by  the 
upper  branch  of  the  parabola  and  its  extension  into  the  imagi¬ 
nary  sound  speed  region.  The  net  of  characteristics  is 
illustrated  by  Fig.  4.  The  sonic  parabola,  which  is  also 
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t 

t 
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shown,  has  in  Lagrangian  coordinates  the  equation 


r-/  /  do3-  r XzL  fs*-s)  (29) 

t  “  srf^jic  V  *  (  y 

For  rt  ference  the  position  in  tine  of  the  x  -  o  Eulerian 
coordinate  has  been  included  in  Fig.  4.  Its  equation  is 

t  -  -%■  s  oo) 


For  a  constant  acceleration  gun  with  the  projectile  at 
s  -  o  for  t  -  o  we  find  from  equation  28  that  at  the  time 


(31) 


sonic  conditions  are  reached  at  the  base  of  the  projectile. 
Sonic  velocity  occurs  at  the  barrel  entrance  x  -  o  at  the 
intersection  point  of  the  parabola  represented  by  Eqs.  (28) 
and  (29) .  We  obtain 


S  = 


Gp 

3  -yj 


(32) 


t 


/  ± 

o<  /  3-r 


(33) 


The  Lagrangian  presentation  illustrates  best  the  fact 
that  the  forward  and  backward  running  waves  completely  cancel 
each  other  and  that,  consequently,  the  properties  of  the  flow 
cannot  change  with  time. 

The  effort  of  drawing  the  characteristic  net  in  Lagrangian 
coordinates  can  be  reduced  greatly  by  using  an  a  versus  u  plot. 
The  characteristics  in  these  coordinates  are  then  two  families 
of  straight  lines  with  slopes  of  -  •  This  is  in  general 
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true  independent  of  the  particular  flow  under  consideration. 
The  unique  difference  here  is  that  for  a  constant  acceleration 
flow  in  Lagrangian  coordinates  the  sound  speed  is  a  function 
of  the  initial  coordinate  alone.  We  have  then  u  -  u(t)  and 
a  -  a(s)  and  thus  an  a,  u  plot  also  presents  an  s,  t  plot. 

In  nondiaensional  coordinates  the  relations  are  u/c(- 
»nd  «/*. O  ~  ** * 

The  straight  line  characteristic  net  is  shown  in  Fig.  5 
for  the  case  of  .  The  sonic  parabola  degenerates  into 

two  straight  lines  as  shown.  The  Eulerian  coordinates 

x  -  o  prescribe  a  hyperbola  of  the  equation 

(a/ W - V"  ~  /  (34) 

Point  A  gives  the  time  when  at  the  base  of  the  projectile 
and  point  B  when  at  the  barrel  entrance  sonic  conditions 
occur.  The  C+  and  the  C_  characteristics  through  point  B  are 
the  saae  characteristics  which  are  shown  in  Fig.  2.  The  C_ 
characteristic  through  point  B  is  tangent  to  the  hyperbola. 
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